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I. INTRODUCTION 



As the only charged gauge boson in the Standard Model (SM) whose mass is generated 
through the electroweak symmetry spontaneous breaking (EWSB), the precise measure- 
ments on W boson's mass and decay width offer a way to distinguish the different EWSB 
mechanisms, and effectively put a stringent constraints on the Higgs boson mass and New 
Physics beyond the SMpJ. Furthermore, the study of the productions of W bosons can be a 
probe for searching New Physics since the possible new particles may decay to W bosons. By 
reason of that, recently, there have appeared several papers discussing "W boson physics", 
at Tevatron, ILC and LHC, respectively j^j. 

Before the LHC era, only two rare decays: W + — > 7r + 7 and W + — > Dfj, were experi- 
mentally studied, and the respective upper limits for their branching fractions are 8 x 1CT 5 
and 1.3 x 10 -3 {3] 1 . The LHC has been running for two years. Due to its high luminosity, the 
huge number of W bosons are and will be produced. According to the estimations in jl, 5, 6|, 
the production cross section of W boson is around ~ 10 2 nb at the LHC with the center en- 
ergy 7 TeV or 14 TeV. It is expected that hundreds of million W bosons will be produced per 
year at the LHC with the conservative expectation of the luminosity ~ 10 32 cm _2 s _1 (which 
is twp orders of magnitude smaller than the desired LHC luminosity ~ 10 34 cm _2 s~ 1 ). This 
makes the LHC also a VT-boson factory which facilitates the precise experimental study on 
W boson physics, especially W boson rare decays. 

In this article, we will discuss a class of W boson inclusive decays to the S-wave quarko- 
nium and B c meson, such as rj c , J/ip, r/b, T, B c and B*. Because W boson is so heavy 
and participates merely the electroweak interactions, the study on W boson inclusive de- 
cays to quarkonia and B^ mesons 2 , offers another great place for the precise test of the 
perturbative quantum chromodynamics (PQCD) and the quarkonium production mecha- 
nism. As a quarkonium is presumed to be a non-relativistic bound state of heavy quark and 
anti-quark, a wonderful theoretical tool to deal with the processes involving quarkonium is 
non-relativistic quantum chromodynamics (NRQCD), in which the low-energy interactions 



1 Recently, the CDF collaboration just set a new 95% confidence level upper limit on the relative branching 
fraction r(TT + — > /T(W + — > e + v) at G.4 x 10 which is a factor of 10 improvement over the previous 
limit 



2 As in §[ 0, Q|, we treat as a non-relativistic bound state of c and b in this paper. 
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are organized by the expansion in v, the typical relative velocity of the heavy quark and anti- 
quark inside of quarkonium pj . The inclusive production cross sections can be written as the 
product of the perturbatively calculable short-distance coefficients and the non-perturbative 
NRQCD matrix elements. At the leading order in v, the only non-perturbative NRQCD 
matrix element of the S-wave quarkonium is proportional to the Schrodinger wave function 
at the origin squared. 

On the other hand, since the W boson mass m w is much greater than the heavy quark 
masses TOf, iC , the parton fragmentation may dominate the W boson inclusive decays to 
quarkonia and B^ mesons. The pioneer works in this field were done two decades ago, 



the universa 



m 



0, Q, 



fragmentation functions for different quarkonium and B c meson were given 



By these universal fragmentation functions, the direct production rate of 
quarkonium and B c meson at large transverse momentum in any high-energy process can be 
approximated to be the product of the parton-level production rate and the universal frag- 
mentation probability, and tbe corresponds computation is m uch easier than that done by 



within the NRQCD framework. In Ref.[8 



10[ , for the Z° inclusive decays to the S-wave char- 



monium, the authors did find great agreement between the results obtained by the parton 
fragmentation approximation and the complete leading order PQCD calculation. 

However, in some processes, the fragmentation may not be the leading process, but main 
process when a certain condition of the fragmentation mechanism is not sufficiently fulfilled. 



In Ref. 



12| it was found that the fragmentation contribution is important but not dominant 



for top quark decays into quarkonium. Thus, it is worth examining if the fragmentation 
mechanism works in W boson inclusive decays to quarkonium and B c meson. 

In the following sections, we present the complete leading order calculation of W + bo- 
son decays to quarkonium r] C} J/ip, rjf,, T and B^ + mesons; then we compare our results 
with that obtained by the fragmentation approximation; finally we discuss the theoretical 
uncertainties in our calculation. 



II. FORMALISM 

Some typical Feynman diagrams for quarkonia and B^ mesons hadronic production 
through W + decay are shown in Fig. [TJ We will calculate them at the leading order in a s 



3 



Vc(or J/0) 



rp,{ar T) 




(*) 

FIG. 1: The typical Feynman diagrams for quarkonia and Be productions from W decays. Some 
diagrams with the different gluon attachment are not shown. 

and v. In NRQCD, a non-relativistic bound state of heavy quarks Q and Q' is considered 
as an expansion of a series Fock states. The leading Fock state for the S-wave hadron is 
constructed as 



H[ 2s+1 Sj}(p,X)) 



where denotes the mass of hadron and all the states are relativistically normalized. Here 
p is momentum of hadron, i,j stand for the color indices which run from i,j = 1, N c with 
N c = 3 for QCD, and C(J, A; 1/2, Ai, 1/2, A2) is the C-G coefficient with A, Ai and A2 being 
the third components of spin indices. The non-perturbative parameters ^2s+i Sj (pq) is the 
Schrodinger wave functions in momentum space. Thus, to project out the amplitude for 
the S-wave hadronic states from the complete patron level amplitude at the leading order 
of v, practically we do the following replacements for the heavy quark Q and -anti-quark Q' 
spinors 

•*«'>^> - + m ° + mQl) 7W: (2) 
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Here e M is the polarization vector of 3 Si state, ^>2 S +i Sj (0) is the Schrodinger wave function at 
the origin. At the heavy quark limit, we set \l/3 5l (0) = \l/is o (0). Throughout the paper, we 
adopt the relation van = itlq + itlqi since we are doing the LO calculation in v. Obviously, 
we have rriQ = vtiqi for quarkonia, and itlq = mb, tuqi = m c for B^*\ 



A. W + ->■ rj c (or J/ijj) + c + s 



There are two Feynman diagrams for W + boson decay into a S-wave charmonium state 
associated with c and s. Implementing the Feynman rules of the Standard Model and 
projectors in Eq.©, the amplitude for W + — » n c + c + s is 

167T qa s V cs ^^(O) . . . . 



where 



A 1 = u(p 6 ) 7 a 7 5 (2m c + ^(l- 7 
(-fa -fa-fi 



l a v{ph)i ; tz, (5) 



(m c + fa + 
(P3+P4+Pe) 2 ~ ml' 



A — / \ a 5/q . in (m c + fa+fa+fa) / 



with p, and ^6 being the momenta of n c meson, s quark and c quark, respectively; 
and the momenta of c quark and c quark in n c meson; V cs the CKM matrix element; 
g = ej sin^vi/ with 9w the Weinberg angle and e the unit electro-charge; a = e 2 /(4n) and 
a s = g 2 /(4ir) are the electromagnetic and the strong coupling constants, respectively; ew 
the polarization vector of W + boson. In this paper, we set u, d and s quarks massless, and 
ignore the dependence on the relative momentum between quark and anti-quark in meson, 
and set P3 = ^4 = p/2. The spin-averaged partial decay width of W + —> r\ c + c + s reads 

dT = ^3 E W^d** (7) 

spins 

where m w is the mass of W boson, s\ = (p + p^) 2 and s 2 = (jps +Po) 2 - The explicit analytic 
expressions for the square of the amplitude for W + — > rj c + c + s is given in Appendix A. 

For the complete leading order calculation of W + — > J ftp + c + s, we just repeat the 
calculation that had been done in Ref.[9|. We find that there is a misprint in the term 



n 

ixy 3 r'My^(l + r — r' — r)~ 2 in I ^2 1 2 given in Appendix B of [9], which should be 



corrected to —8xy r'M^(l + r — r' — r) . 
B. W + -> rj b (oi T) + X 

W + — > ^(ot T) + c + b is the leading order process of bottomonium production through 
W + decay in the expansion of a s . However, such processes are the CKM-suppressed with the 
corresponding CKM factor V c b = AX 2 in the Wolfenstein parameterization where A = 0.813 
and A = 0.225. Numerically, the Wolfenstein parameter A ~ a s (2m c ) = 0.26. Hence, we will 
consider some processes that are higher order in expansion of a s but lower order in expansion 
of the Wolfenstein parameter A. For instance, the process W + —trjb + c + s + g depicted by 
the last diagram in Fig. [TJ is associated with the factor glV U( i ~ g 3 s A which is numerically 
comparable to the factor g 2 V c b ~ g 2 X 2 accompanied with the process W + — > rjb + c + b. 
Similarly, we also consider W + — > T + c + s + 7 for W decays to T. Its amplitude is 
proportional to V^e 2 , which is numerically comparable to the factor g 2 s V c b accompanied 
with W + ->■ T + c + b. 

The calculation of rjb production is similar to the r] c production described in the previous 
subsection. Here we just show the results for the production of T. The complete amplitude 
for W + decay to T with c and b is 

167T qa s V c b ^t(O) , . . s ,„ s 



with 



Ai = u{p 5 7 - — — -/ w (l-7 ^ T 

(P3 +P5 +P6) ~ m c 



(2m b + fylavipe)-, — vj, (9) 

{P3+P6) 

A 2 = U{p 5 )^ W (l-J 2^ FT 

[P3 + Pi + Pe) ~ ™ h 

(2m b + 0h a v(pe)- ■ (10) 

{P3+P6) 

where p, p§ and p§ are the momenta of T, c quark and b quark, respectively; ^3 and p^ are 
the momenta of b quark and b quark in T meson; ex is the polarization vector of T boson. 
The explicit expression for the square of amplitude is given in Appendix B. 
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FIG. 2: Relevant electroweak Feynman diagrams which contribute in the production of T. 



Because of the CKM suppression, the electroweak contributions depicted in Fig. [2] is 
required for the production of T as we argued above 3 . Taking W + — > T + u + d for example, 
the amplitude for this process is written as 

M = V3ge 2 V ud -^A(A 1 + A 2 + A 3 ), (11) 



with 



At 



A, 



iy[ 7Q &(2m 6 +j*)] 

9v^(p 3 +P4) 2 

x^ 6 ) 7 " / 3 +^ 6 Ml - 7 5 )^(p 5 ) , (12) 

(P3 +Pi + P&) 



18^3 + Pi ) 2 



A, 



xu(pe)Ml ~ h M \ Q -fvfa) , (13) 

{P3+P4+P5V 

Tr[ 7 ^(2m 6 + $]u(p 6 )y(l - 7 5 )e(>(p 5 ) 



X ((Pfi + P6 ~ P)n9au ~ (2p 5 + 2p 6 + p) a ^x/ 
+ (Pfi + P6 + 2p) (14) 



where p, ps and p6 are the momenta of T, d quark and w quark, respectively; p 3 and ^4 are 
the momenta of b quark and b quark in T meson. Here again we set p 3 = p^ = p/2 in the 
calculation. 



3 However, for that of J ftp, the electroweak contribution can be ignored in LO, which is two orders of 
magnitude less than its QCD contribution. 
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W + — > r]b + c + s + g dominates the decay W + — >■ r] b + X. Its four-body decay am- 
plitude is quite complicated, so we do not present the complete analytic expressions here 
but the numerical results in the next section. The potential infrared-divergences arising 
from the regions where the momentum of the gluon becomes soft vanish due to the "color- 
transparency" , and the potential collinear-divergences arising from the region where the 
momentum of the gluon becomes collinear to the r\ b and c-quark (or s quark) are regulated 
by the b quark mass. Thus, the decay width of W + — >■ rj b + c + s + g is infrared-safe. 

C. W + -»• B { * )+ + b + s and W+ -»• B { * )+ + c + c 

There are two processes, namely W + — > + b + s and W + — > + c + c, for 

be bound states production in W + decay at the leading order. W + — > + g + g and 

w + B { ; )+ + 7 can also contribute, however their contributions vanish at the large m w 
limit, and thus are almost two and three orders of magnitude lower than W + — > + b + s 
and W + — > B^ + + c + c numerically. 

For W + — > + b+s, there are two Feynman diagrams and the corresponding amplitude 

is 

M = ^sVcs {a a) (15) 

3 2v/2 2^3(m c + m b y h K J 

with 

[P3 + Pa) 

-7^(P5j, (16) 



(pa + P5 + Pe) 



2 



M (Pe)7 a 7 5 ( m c + m 6 + ^)7 a 



(P3 +Pe) 2 
(m c + ^3 + ^4 + i>% 

* 2 c 



Ml-T 5 )^)- (17) 



(P3 +P4 + Pe) 2 - ™c 

Here p, ^3, P4, P5 and p$ are, respectively, the momenta of B£, b quark and c quark in 
B+ meson, s quark and b quark. The corresponding amplitude at the leading order for 
W + ->■ 5+ + c + c is 

M = ^^L(A 1 + M (18) 
3 2v/2 2v/3(m c + m 6 ) V 1 v ; 
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with 



A = u(pe)7 7 ; ; ^ j iv 1 - 7 

l b (m c + m h + f}^ a v(p b )- ■ — , (19) 

KPa+PsY 

a 2 = «^M-T7-rri2 2-7 

7 5 (m c + m b + ^)7 Q w(p 5 )- ■ — , (20) 

(Pi + PsY 

where p, p 3 , p 4 , p 5 and p 6 are, respectively, the momenta of B c , b quark and c quark in B c 
meson, c quark and c quark. 

Similarly, we can obtain the amplitudes for B* production. In Appendix C, we detail 
the square of the amplitudes for both W + —y Be + b + s and W + — > B c + c + c. 



III. NUMERICAL RESULTS AND ANALYSIS 
A. Input parameters 

Now we can employ the above formalism to evaluate the decay widths of W + boson to 

(*)+ 

quarkonia and Be mesons. We adopt the mass parameters for c quark, b quark and W 
boson as follows 

m c = 1.50 GeV , m b = 4.90 GeV and m w = 80.4 GeV. 

The Schrodinger wave functions at the origin for J/ip and T are determined through their 
leptonic decay widths T ee , at the leading order in both a s and v, which is Q, 

kt/ rmi 2 - m2 ^) Tee 

167ra 2 e2 (b) 

Taking Vf e = 5.55 keV, Tj e = 1.34 keV, a = 1/137 and e c = 2/3, we obtain 
|^(0) | 2 = 0.0447 GeV 3 and |^ T (0)| 2 = 0.403 GeV 3 . 

In principle, we can extract ^^(O) from the decay width of r] c — y 77. However, since 



experimental data uncertainty is still large, which is B{j] c — y 77) = (6.3 ± 2.9) x 10 



-5 



he 
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TABLE I: Decay widths and branching fractions of quarkonium and B c meson through W + inclu- 
sive decays. Here the qiijj represents one kind of cs, bs, cc, cb, ud, H + Vfi and csg; the QiQj one 
kind of cs and cb. 



W+^Hq^ T{W + -> i? % ^)(keV) f/(W+ -> ff ft ^)/T(W r+ -> QiQ;) 



iy+ - 


> J/tpcs 


90.4 


1.27 x 10~ 4 




-> r/ c cs 


87.5 


1.23 x 10~ 4 


- 


-»■ B+6s 


6.3 


8.84 x 10~ 6 


- 


> S*+6s 


5.4 


7.61 x 10~ 6 


w + - 


-> B+cc 


0.4 


3.64 x 10~ 4 


w + - 


> B*+cc 


0.5 


5.23 x 10~ 4 


w + - 


-> 


0.015 


1.41 x 10~ 5 


w + 


-> Tc6 


0.016 


1.49 x 10~ 5 


w + - 


-> Tud 


0.014 




W + —7 


- T>+z^ 


0.013 




w + - 


-)• Tftcsg 


0.517 





the corresponding extraction of ^^(0) is not applicable in our numerical analysis. Instead, 
we set ^^(O) = ^jm(O), which is the consequence of the heavy quark spin symmetry in 



NRQCD at leading order v\L\. Similarly, we adopt * % (0) = ^ T (0). 



For Be , we apply |Wr 



Buchmiiller-Tye potential 
adopted, i.e. 



0.1307(GeV) 3 given in Q by using the 



15] . In addition, for consistency the leading order a s running is 

47T 



;n-2V3)in(u 2 /A^)' 

where we take Aq CD = 200MeV and rif = 4 as in 14J. We choose the typical renormaliza- 
tion scale fi = 2m c for charmonium and production, and fi = 2m& for bottomonium 
production, correspondingly a s (2m c ) = 0.26 and a s (2m&) = 0.19. 
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B. Decay widths and feasibility at LHC 



By the input parameters given above, we list the partial widths of W decays to the 
S-wave quarkonia and mesons in Table HI Employing these partial widths, one can 
estimate the event numbers of quarkonium production through W decays at the LHC. 
Considering that the LHC runs at the center-of mass energy = 14TeV with the luminosity 
£p_ p = 10 34 cm _2 s _1 , and the cross section aw+ = 10 2 nb at the LHC {f5|, the number of W + 
events per year is expected to be 3.07 x 10 10 . Then we present the event rates of W decays 
to quarkonia or Be in Table [III 

The most readily identifiable quarkonia are J/ip and T, because their leptonic decays 
have clear signals and relatively large branching fractions. With B(ip — > fifi) = 5.93%, 
B(T — > /U/i) = 2.48% and B(T — > rr) = 2.60%, we list the expected di-leptonic signals from 
J/ijj and T decays in Table II III For experimental detection for such decays at the LHC, one 
has to look for the relevant events pp — > W + X with W — > J /ip([i + ^) + 2 jets. According 
to [if]], the efficiency of reconstruction of J/ip from its dimuon decay channel is around 40% 
for pp collision at LHC experiments. However, to tag whether such J /if) event is really from 
W + decay, one must reconstruct W + from + 2 jets event. Therefore, the additional 

information about relevant QCD background of such events could be crucial as well. 

The most promising decay channels for reconstruction rj c are r] c — > KKtt , rjiTTc and etc, 
which branching ratios are around a few percent. Thus, it should be possible to observe 
W + decays to r] c at the LHC. Certainly, such measurements also suffers the difficulty of 
reconstruction of W + as mentioned above. 

B£ — > J/ip7T + is the most promising channels to identify B+ at colliders. However, 
considering their small branching fractions and the efficiency of the event reconstruction, 
the measurements on W decays to B c at the LHC would be quite difficult. With the similar 
reasoning, it would be barely possible to measure W decays to r/b at the LHC due to the 
lower branching fraction and the experimental difficulty for identifying rjt,. 
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TABLE II: The expected annual event numbers of the quarkonium and B c meson from W + decay 
at the LHC with yfi = 14TeV and C p - P = 10 34 cm- 2 s _1 . 







Vc 


B+ B*+ 


'lb 


T 


Events number(xl0 4 ) 


133 


129 


9.8 8.8 


1.5 


0.12 


TABLE 


III: The corresponding annual leptonic events of 


J/ip and T. 












T 




Decay channel 


e + e~ 




e + e 




r + r - 


Events 


7.9 xlO 4 


7.9xl0 4 


30 


30 


31 



C. Comparisons with the fragmentation mechanism 

Since m w ^> , m c , one can argue that the fragmentation mechanism may be the 
dominant contribution to the inclusive decay rate of the W + into the quarkonium and B c 
meson which survives in the limit m w /mb tC — » oo. In the fragmentation mechanism, the 
hadron H with energy E is produced by the fragment of a type i parton with energy E/z 
(z is the longitudinal momentum fraction of H relative to type i parton) directly from the 
decay of W + boson. The possibility of % — > H is presumed to be described by the universal 
fragmentation function D^niz)- Putting all the possible parton fragmentation together, 
the differential decay width of H inclusive production can be written as 

dT(W + -> H(E)+X) 




where dT is the parton level differential decay width, and the sum is over the parton type i 
and the hadron H's longitudinal momentum fraction z. At the leading order of a s , dT does 
not depend on the longitudinal momentum fraction z. Thus, the total decay width turns to 
be 

T(W + -)■ H{E) +X) = J2 -)• i(E/z) + X)P^ H , 

i 
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Pi->H = 



dzDj 



(22) 
(23) 



where Pj-^ij is the so-called fragmentation possibility. 

After paying the price of 0(rriH/mw) power corrections, Eq. (l2~Tj) has a number of advan- 
tages in calculation. The parton level differential decay width df is easy to be calculated 
and dependent only on the typical short-distance scale m w . The fragmentation functions 
are universal and dependent on the hadronization scale, and therefore, they can be either 
extracted from the experimental measurements or calculated from certain phenomenological 
models. Fortunately, the fragmentation functions for the S-wave quarkonium and B c meson 



can be calculated perturbatively 
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111 ], together with all the non-perturbative hadroniza- 



tion effects being parameterized into the Schrodinger wave functions of the mesons at origin 
at the non-relativistic limit. 



In calculation we take the following fragmentation probabilities from 

32a 2 (2m c )|^,(0)| 2 /1189 
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P( C ->• T] c ) 

P(b B+) 



27m 3 

32« 2 (2m c )|^ 
27m? 



n 30 

2 / 773 
~30 



- 571n2 



37 In 2 



k* 2 (2m c )|^ 



27m 3 
la 2 s (2m c )\^ BS (0)\< 



-f 



m r 



27ml 
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P{b Bl + ) = 

where the function f(r) and g(r) are 

8 + 13r + 228r 2 - 212r 3 + 53r 4 



mi, + m c 
mt, + m c 



(24) 

(25) 
(26) 
(27) 



f(r) 



15(1 -r) 5 

r(l + 8r + r 2 — 6r 3 + 2r 4 ) 
+ n ln ( r ) 



(28) 



9{r) 



24 + 109r - 126r 2 + 174r 3 + 89r 4 
15(1 -r) 5 
r(7 - 4r + 3r 2 + 10r 3 + 2r 4 ^ 



ln(r). 



(29) 



(l_ r )6 

Pb->T an d Pb^-q b can be obtained from (23) and (24) by substituting the mass m^ for m 
And P „(*)+ can be also got from (25) and (26) by interchanging m^ and m c . 
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TABLE IV: The branching fractions for the quarkonium and B c meson production predicted by 
the complete leading order calculation and fragmentation approximation and their comparisons. 









LO calculation (B) 


fragmentation (B*) 


B*—B 
B 


F(W + - 


> J/^cs)/T(W + - 


-> cs) 


1.27 x 


1.31 x 1(T 4 


3.3% 


T{W + - 


-> T] c cs)/r(w + - 


> cs) 


1.23 x 10" 4 


1.27 x 1(T 4 


3.6% 


T(W + - 


» B+bs)/T{W + - 


-)■ cs) 


8.84 x 10" 6 


1.03 x 10" 5 


17% 


T(W + —, 


> B*+bs)/r(w + - 


cs) 


7.61 x 1(T 6 


8.92 x 1(T 6 


17% 


F(W+ - 


+ B+cc)/T(W+ - 


-> cb) 


3.64 x 10" 4 


3.85 x 10- 4 


6% 


F(W+ - 


> B* c + cc)/T(W + - 


-»• cfe) 


5.23 x 1(T 4 


5.41 x 1(T 4 


3% 


T(W + - 


-> rj b cb)/T(W + - 


>■ eft) 


1.41 x 10" 5 


1.76 x 1(T 5 


24% 


T{W + - 


-> Tcb)/T(W+ - 


► c5) 


1.49 x 1(T 5 


1.81 x 10~ 5 


21% 



The numerical results on the decay widths of W inclusive decays to quarkonium and 
B c meson by using the fragmentation approximation are listed in Table IIVI We also list 
the comparisons between the predictions given by the fragmentation approximation and 
the complete leading order calculations. We find that the difference between the predic- 
tions by the two methods is negligible for the charmonium productions. However, for the 
bottomonium and productions, the differences is sizable but understandable, since the 
power corrections to the fragmentation mechanism is order of 2mb/m w , which is numerically 
around 10%. 

D. Theoretical uncertainties 

At last, we investigate the theoretical uncertainties in our calculations for the J/ip and r] c 
production from W boson decay. The main sources of uncertainties include the wave function 
at the origin, the strong coupling constant and the quark mass. After some numerical check, 
we find that the biggest uncertainty in our results is from the choice of the charm quark 
mass. Here we present the dependence of our results on the charm quark mass in Fig. [31 It 
shows that we can get a production rate of charmonium which is 1.6 times as much as we had 
calculated if we adopt m c = 1.27 GeV. Besides, the production rate from fragmentation is a 
slightly bigger than that from complete leading order computation. The difference between 
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the LO calculation and the fragmentation approximation is insensitive to the charm mass. 
The second largest uncertainty is due to the choice of renormalization scale as mentioned 
in the front of this section. Also we lay out their relations in Fig. [3j The dependence 
on the renormalization scale can be reduced by considering the NLO QCD corrections and 
employing the Altarelli-Parisi equations to resum the large logarithms. 




m (GeV) fi (GeV) 

FIG. 3: The partial width of the process r(VF + —> J/ipcs) versus the mass of charm quark(upper, 
fx = 3 GeV) and the renormalization scale(lower, m c = 1.5 GeV). The dashed line corresponds 
to the prediction by the complete LO calculation, while the solid line corresponds to that by the 
fragmentation approximation. 



IV. CONCLUSIONS 



In this paper, we show that there will be around 10 6 J/^S and rj c , 10 4 Be mesons and 10 3 
T and rjb annually produced from W boson inclusive decay at the LHC, if the machine works 
at the center energy y/s = 14TeV and with the luminosity £ p _ p = 10 34 cm _2 s _1 . And the 

fragmentation mechanism works well for W boson inclusive decays to the S-wave charmo- 

(*) 

nium, bottomonium and Be mesons. The numerical calculation shows that the difference 
between fragmentation approximation and complete LO calculation for W boson decays to 
charmonium is around 3%, and it is not magnified by considering the uncertainties of theo- 
retical parameters. But for W decays to B^ mesons and bottomonium, the differences rise, 
but it is understandable. 
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Appendix A: The square of the amplitudes for W + — > rj c + c + s 

Here we introduce parameters 

sio =Pi-P, si 5 =pi-p 5 
with pi, p and p$ being the momenta of W + , i] c and s quark, respectively. 



E|A| 2 

1 

(ml(m 2 c + m 2 w - s 10 ) 2 (m 2 c - m 2 + 2si 5 ) 2 ) 
x (-64(12m^m 2 - m s w + 4s 2 si 5 (sio + s 15 ) 

+4m^(2s 10 + si 5 ) + 2m 2 u s 10 (7s 2 + 4s 2 5 + 10s 10 s 15 ) 

-mt(19s 2 + 20s 10 si 5 + 4s 2 5 ) 

+m*(-5m* - 3s 2 + 24si si5 + 16s 2 5 

+m 2 y (22s 10 + 40s 15 )) - 2m 2 (3m^ - m 4 w (17s w + 14s 15 ) 

+sio(-s? + 10si si5 + 8s 2 5 ) 
+m 2 y (21s 2 + 28s 10 si5 + 8s 2 5 )))), 

S|A 2 | 2 

1 

{m 2 w {m 2 c - m 2 + 2s 15 ) 4 ) 

x(64(27m 6 c m 2 w + (m 2 - 2s 15 ) 2 (mt - 2m 2 (s 10 - 2s 15 ) 
-4s 15 (si + s 15 )) - m 2 (m 2 - 2s lb ){7m 4 w 

-8si 5 (5Si + S15) - 2TO 2 (l0si + S15)) 

-3m^(7mt + 4(3s 10 - 5si 5 )si 5 + 2m 2 (3s 10 + 8si 5 )))) , 
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ERe(AiA 2 *) 

1 

(ml(m 2 c +m 2 w - s w )(m 2 c -m 2 w + 2s 15 ) 3 ) 

x (64(9m 6 c m 2 w + (m 2 w - 2s 15 )(2m 6 w + 4si si 5 (sio + s 15 ) 

+2m 2 w s 10 (3s 10 + 2s 15 ) - ml,(7s 10 + 4s 15 )) 

+m 4 c (2s 15 (21s 10 + 4s 15 ) + 5m 2 w (3s 10 + 8s 15 )) 

+m^(5m^ - 16m^si 5 + m^(-14s^ 

v + 44) + 4s 15 (-54 + 4s 10 s 15 + 44,)))) . 

Appendix B: The square of the amplitudes for W + — > T + c + b 

Here we define parameters: 

SlO =Pl-P, S15 = pi • J9 5 

with pi, p and p 5 being the momenta of W + , T and c quark, respectively. 

f 

(m^m 2 -m 2 c + m 2 w - S\^ 2 {-m\ + m 2 c + m 2 w - 2s i5 ) 2 ) 
x (64(66m^m 2 — 2m\m 2 w + — 8m^,si + S5m w sl 

-46m 2 w sl - 4m£,si5 + 20mtsi si 5 - 20m 2 u sl s 15 

-4s? si5 + 4m* s? 5 - 8m^ 10 s^ 5 - 444 

-mc( 3m t + m l(22si - 8si 5 ) - s w (7s w + Ss 15 )) 

+ml(A2m 2 m 2 w + m A w + 34 - 32si si5 + 164 

-2m 2 w (97s l0 + 32s 15 )) - 2m 2 (m 4 w (15s 10 

-2s 15 ) + (-9s? - 28s 10 si 5 + 44) 

+sio(lls? + 18s w s 15 + 8s? 5 )) - 2m 2 b (5m 4 c m 2 w + 2m° 

+4( s io - 10si 5 ) + mt(27si + 10si 5 ) 

-m 2 w (854 + 36s 10 s 15 + 44) - m;;(25mt + 234 

+28si si5 + 84 - Am 2 w (7s 10 + 5s 15 ))))), 
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S|A 2 | 2 

1 

(m 2 w (-m 2 + m 2 c + m 2 w - 2si 5 ) 4 ) 
x(64(27m£m 2 - (m 2 c + m 2 w - 2s 15 ) 2 {m 2 c m 2 w 

-m* + 2m 2 w (s 10 - 2s 15 ) + 4si 5 (si + s 15 )) 

-ml(m 2 c + m 2 w - 2si b )(3m 2 c m 2 w + 7m* - 8si 5 (5si + s 15 ) 

-2m 2 (10s 10 + Sl5 )) + m 4 (73m 2 m 2 - 21m* 

-6m 2 w {3s w + 40s 15 ) + 4s 15 (-9s 10 + 47s 15 )))), 

1 

(m^m 2 , -m 2 c + m 2 w - s 10 )(-m 2 + m 2 c + m 2 - 2si 5 ) 3 ) 
x (-64(-51m^m^ + m 4 (-49m 2 m 2 + 40m* 

+2(13s 10 - 44s 15 )s 15 + 7m 2 y (13s 10 + 16s 15 )) 

+m 2 b (3m 4 c m 2 w + 25m 6 w - 4m 4 y (7s 10 + 6s 15 ) 

-2m 2 (2mt + 9sj + 22m 2 s 15 - 2s 10 s 15 - I2s 2 15 ) 

+m 2 w (-38s 2 10 - 24s 10 si 5 + 4s 2 5 ) 

+4s 15 (-s 2 + 12s 10 si 5 + 4s 2 5 )) 

+ (m 2 +m 2 w - 2si 5 )(m 4 m 2 y + 2m 6 w + 4s w s 15 (s w + s 15 ) 
+2m 2 w Sio(3si + 2s 15 ) - m 4 (7si + 4s 15 ) 
+m 2 (3m 4 + 2s 10 (s 10 + s 15 ) - m 2 w (3s w + 2s 15 ))))). 

Appendix C: The square of the amplitudes for Be production 

First we present the squared matrix elements for W + —¥ B+ + b + s. We introduce 

SlO =Pl-P, Si 5 =Pl-p 5 

with pi, p and p^ being the momenta of W + , and s quark, respectively. 
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1 



m\m 2 w (m h (m 2 c + m 2 w ) + m c (m 2 c + m 2 w - 2s w )) 2 

X (ml -ml + 2s 15 ) 2 (16(mb + m c) 2 {2mtm c ml{ml 

-m 2 w - 2s 15 ) + ml(m 4 m 2 w + m 2 w (-m 4 w + 2m 2 w (s w - 2s 15 ) 

+8si 5 (sio + si 5 )) - 2m^(2si 5 (si + S15) 

+3m 2 w (3s 10 + 4si 5 ))) - m^(m^ + - 2si )(3m^m^, 

-ra£, + 4m^(si + S15) - 8si si 5 (sio + s 15 ) 

-Am 2 w (s 2 10 + S10S15 + s? 5 ) + m l(- 2m i + 4m l(sio + si 5 ) 

+4(-Si + S10S15 + Si 5 ))) - 2m h mJ(lm\m 2 w 

+mt(-m 4 w - As 2 10 + 12s 10 si 5 + 8s 2 15 + 2m 2 w {s w + 6s 15 )) 

-m 2 w (m 2 w - 2s l0 )(m1 - 2m 2 w (3s 10 + 2s 15 ) 

+4(2si + 3s 10 s 15 + «? 5 )) - m 2 c (5m 6 w + 8si Si 5 (3si + 2s 15 ) 

-2m* (llsi + 8si 5 ) + 4m 2 w (6s 2 10 + 6s 10 s 15 + s 2 15 ))) 

+m\(-22m\m 2 w + m\{\\m 4 w - 8s 15 (5s 10 + 3s 15 ) 

-4m 2 w (7s w + 13si 5 )) + 2m 2 c (5m 6 w + 8si si 5 (sio + s 15 ) 

-mt(29s 10 + 22s 15 ) + 12m 2 w (3s 2 10 + 4s 10 s 15 + s 2 15 )) 

+m 2 w (m e w - 16sio(si + s 15 ) 2 - 2m^,(5si + 2s 15 ) 

+4m 2 w (Qs 2 10 + 7s l0 s l5 + s 2 15 ))) 

-2m\m c (§m\m 2 w + m 2 c (-hm A w + 4s 15 (3s w + 2s 15 ) 

+m 2 w (20s w + 26si 5 )) - m 2 w (m 4 w - I2m 2 w (s w + s 15 ) 

+4(4s? + 8s 10 si5 + 3s? 5 ))))), 



S|A 2 | 2 

1 

(m 2 b m 2 w (m 2 c - m 2 , + 2si 5 ) 4 ) 

x (16(m 6 + m c ) 2 (2mlm c (m 2 m 2 w - m 4 w 

-10m 2 w s 15 + 8s 2 15 ) + 2m b m c (Qm 4 m 2 w - (m 2 w - 2s 15 )(m 4 v 

-8s 10 s 15 - 2m 2 w (2s 10 + s 15 )) - m^(5m* + 8(si - Si 5 )s 15 

+2m^(2si + si 5 ))) + (m 2 - m 2 + 2s 15 )(4m^4 
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-(ml, - 2s V r ) )(m A w - 2m 2 w (s w - 2s 15 ) - 4:S 15 (s 10 + s 15 )) 
+m 2 c (3m 4 w - 4si 5 (4s 10 + 3s 15 ) + m 2 w (-8s 10 + 12si 5 ))) 
+ml(9m 4 m 2 w - (m 2 w - 2si 5 )(m^ - A(s w - 2s 15 )s 15 
-2m 2 w (s w + 7s 15 )) - 2m^(4mt + 2(s 10 - 10si 5 )si 5 
+m 2 w (s 10 + 22s 15 ))))), 

HRe(A 1 A 2 *) 

1 

m\m\(m h (m 2 c + m^) + m c (m 2 c + - 2si )) 
x (m 2 _ m 2 + 2g 15 )3 (~ 16 K + ^c) 2 (2m^m c (m^ 
-3m^ - 2m^ 15 + 2s? 5 ) - m 3 b (3m 4 c m 2 w 
+2m 2 c (3ml, + 2(s w - 2s 15 )s 15 + 5m^(si + 3si 5 )) 
+m 2 w (3ml, + 4si 5 (si + 2si 5 ) - 2m 2 u (3s w + 5si 5 ))) 
-Am 2 b m c (3m A c m 2 w + 2m£, + 2si 5 (2si + s 15 ) 
-2mi(s w + 3s 15 ) + m 2 w (-2s 2 w + 3 S ? 5 ) 
+m*(-2m* + s 15 (7s 10 + s 15 ) + (4s 10 + llsis))) 
+m c (-((m* - 2s w )(m A w - As 2 15 )(m 2 w - 2(si + s 15 ))) 
+m 4 c (ml, - 2m 2 w (s w - s 15 ) - 4si 5 (3si + 2s 15 )) 
-2m^(m^(s 10 - 2s 15 ) + As 15 (-3s 2 10 - s 10 s 15 + s 2 15 ) 
+m 2 w (-2s 2 w + 8s w s 15 + 2si 5 ))) + m b (-5mlm 2 w 
+m 4 c (3m 4 w - 2m 2 w (s 10 + 2s 15 ) - 8s 15 (5s 10 + 2s 15 )) 
-m 2 w (m 2 w - 2s 15 )(3m^ - 2m 2 w (hs w + 4s 15 ) 
+4(2si + 3s w s 15 + Si 5 )) + m 2 c (m 6 w - 6m^,(2si + s 15 ) 
-Ss 15 (-2s 2 10 + 3s 10 s 15 + 2s 2 15 ) + Am 2 w (As 2 10 
+5s l0 s 15 + 4s 2 15 ))))). 

For W + -> 5* + + 6 + s, we introduce 

SlO =P1-P, S15 =Pl-p 5 



20 



with pi, p and p 5 being the momenta of W + , B* + and s quark, respectively. 

1 

mlml,(mb(ml + m^) + m c (m 2 c + m^ - 2si )) 2 
x (rn 2_ rn l Tol A 2 ( 16 ( m ^ + ^c) 2 (6m^m c m 2 (m 2 - m 2 
-2si 5 ) + 2mlm c (A0m 4 c m 2 w - m^(3mt - 48s 2 - 48si si 5 
+4s 2 5 + 4m 2 y (3s 10 + S15)) - ™ 2 c (5m 4 w + 4(5s 10 - 2s 15 )s 15 
+2m 2 (50si + 23si 5 ))) + ml(?>bm A c m 2 w + m 2 w (-3m^ 
+8(3si - si 5 )si 5 + m 2 (6si + 4si 5 )) - 2m 2 c (8m 4 w 
+2(3s 10 - si 5 )sis + m 2 (27s 10 + 28s 15 ))) + m 2 (m 2 + m 2 
-2si )(9m^m 2 u + m£, - 4m^(si + s i5 ) + 8si Si 5 (sio + s 15 ) 
+4m 2 (5s 2 + S10S15 + «? 5 ) - 2m 2 (mt + 2m 2 (7s K ) + S15) 
-2(s 2 - S10S15 + s 2 5 ))) + 2m b m c (23m 6 c m 2 w + m 2 c (-m 6 w 
+2As 2 l0 s 15 - 2m 4 w (21s 10 + 8s 15 ) + 12m 2 (10s 2 + 2s 10 s 15 
+s 2 15 )) + m 2 w (m 2 w - 2s w )(m 4 w - 2m 2 (3si + 2s 15 ) 
+4(6s 2 + 3s 10 s 15 + s 2 l5 )) + m 4 (9m 4 w - 2m 2 (47s 10 + 6s 15 ) 
+4(s 2 - 3si s 15 + 2s 2 15 ))) + m 2 b (88m b c m 2 w + m 4 c (llm 4 w 
+24s 15 (-2s 10 + si 5 ) - 4m 2 u (72s 10 17si 5 )) + m 2 w (ml 
-2m 4 w (5s w + 2si 5 ) - 16si (3s 2 + 2si si 5 + s 2 5 ) 
+4m 2 (10s 2 + 7si s 15 + s 2 5 )) - 2m 2 c (2m 6 w - 8s 2 l0 s 15 
+mt(43si + 18s 15 ) - 4m 2 (33s 2 + 14si si 5 + 2s 2 5 ))))), 



S|A 2 | 2 

= (m|<(m| -ml + 2s 15 y) + ^nemlm c (m 2 c m 2 w 

-ml - lOmlsm + 8s 2 5 ) + (m\ -m 2 w + 2s 15 )(2m 4 m 2 w 
-(m 2 - 2s 15 )(m^, - 2777^(510 - 2si 5 ) - 4si 5 (si + s 15 )) 
+m 2 c (m 4 + 4s 15 (-2si + s 15 ) - ^m 2 w (s w + 2s 15 ))) 
+2m b m c (4m 4 m 2 u - (to 2 - 2si 5 )(to^ - 8S10S15 - 2m 2 (2s w 
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+s 15 )) - m c 2 (3m* + 8(si - 3si 5 )si 5 + m 2 w (As w + 30si 5 ))) 
+m\{\\m\m 2 w - 3m^, + 6m^si + 20m 2 s 2 5 - 8s 2 5 (3,Sio 
+2s 15 ) - 2m 2 c (Am 4 w + 2(3s w - 22s 15 )s 15 
+m^(3s 10 + 58si 5 ))))), 



ERe{A 1 A 2 *) 

1 

mlm^m^ml + m^) + m c (m 2 c + m^ - 2si )) 

x 7~^ — rr(-16(m fe + m c ) 2 (6mlm c (m 2 c m 2 w - 3m^ 

(m 2 . - ml + 2s 15 ) 3 

-2m^i 5 + 24,) + m 3 b (31m 4 c m 2 w - 2m 2 (25mt + 6(s 10 
-4si 5 )si 5 + 3m 2 (5si + 7s 15 )) + m^(-9m* - 4s 15 (3s 10 
+2s 15 ) + 2m 2 (9si + 7s 15 ))) + Am 2 b m c (10m 4 c m 2 w 
-7m 6 w - 2s 2 15 (2s w + si 5 ) + 2m^(5si + 6si 5 ) 
+m 2 (6s 2 - 4si si5 - 7s 2 5 ) - m^(6m* + (5si - 17si 5 )si 5 
+m 2 (20si + 23si 5 ))) + m c (Am 6 c m 2 w - (m 2 w - 2si )(m^ 
-As\ b )(m 2 w - 2(si + si 5 )) + m^(5m* - 4(si - 2si 5 )si 5 
-14m 2 (si + si 5 )) - 2m 2 c (-6m 2 w (s 10 + si 5 ) 2 
+3m^,(si + 2s 15 ) + 4si 5 (-s 2 + 3si Si 5 + s 2 15 ))) 
+m h (2\m\m 2 w + ml(7m 4 w + 8si 5 (-2si + 5si 5 ) 
-2m 2 (29s 10 + 32s 15 )) - m 2 (m 2 - 2s 15 )(3m* 
-2m 2 (5si + 4s 15 ) + 4(2s 2 + 3si si 5 + 4>)) 
+m 2 (-13m° + mt(4s 10 - 2s 15 ) - 8s 2 5 (7s 10 + 2s 15 ) 
+4m 2 o (10s 2 + 13 Sl0 Si 5 + 4 S 2 5 ))))). 

Then, we present the squared matrix elements for W + — > + c + c. Here we introduce 

sw =Pi-P, si 5 =pi-p 5 
with pi, p and ps being the momenta of W + , and c quark, respectively. 

si^i 2 
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1 

(mjj + m\m c — m\ + m c ml, + m b (-m 2 c + m 2 - 2si )) 2 
1 

m 2 m 2 (m 2 , - m 2 + m 2 - 2s 10 - 2si 5 ) 2 
x(— 16(m b + m c ) 2 (?>m\m 2 w + \^m\m c m 2 w 

-2m\m c (§m\m 2 w - 19777* + 8s 2 - As l0 s l5 - 8s 2 l5 

+Um 2 w (s 10 + 2s 15 )) + 77?f (14m 2 m 2 + 9m* 

-2771^(5510 + 6si 5 ) + 4(-s 2 + s 10 s 15 + s 2 5 )) 

-2mlm c (l§m A c m 2 w - 7m^ + m* (26si + 8si 5 ) 

-4m^(si + 12si si5 - s 2 5 ) + 8sio(-s 2 + S10S15 + 2s 2 5 ) 

-2m 2 c (9m 4 w - 2s 2 + 2m 2 o (s 10 - 10s 15 ) + 4s w s 15 

+8s 2 5 )) - m 2 (30mX, - m 4 (23mt - 12s 2 

+2m 2 (sio - 34s 15 ) + 44s 10 si 5 + 28s 2 5 ) 

+ (777^ - 2s 10 )(m 2 - 2s 15 )(777t + m 2 w (4s 10 - 2s 15 ) 

-4Sio(si + S15)) - 4777 2 (3777^ - 12m* Si 

-4si Si5(sio + s 15 ) + m 2 (3s 2 + 20si Si 5 - 6s 2 5 ))) 

+7T7 2 (-7777^777 2 - 777 2 (m 2 - 2Si 5 )(m* + 8Si Sl 5 
-2777^(810 + S15)) + 777 2 (-777^ + 85^ - 4m* (si - 2Si 5 ) 

+4m 2 (s 2 + 2 Sl0 s 15 - 2s 2 5 )) + m*(9mt 
-2m 2 w (3s 10 + lOsis) + 4(3s 2 + 3s w s 15 + s 2 5 ))) 

-2777 b 777 c (ll7T7^777 2 + TO^(m^ - 2Si )(m 2 - 2si 5 )(m 2 u 
+2Si - 2Sl 5 ) - m 2 (777^ - 8S 2 1Q S 15 + 777* (-10Si + 8S15) 
+4777 2 (2S 2 + 4si Sl5 - 3fi 2 5 )) - 777*(llm* 

-2t77 2 o (3s 10 + 14s 15 ) + 4(s 2 + 5s 10 si 5 + 2s 2 5 ))) 
+m*(-44m*m 2 u + 5777° - 4mt(5s 10 + 2s 15 ) 
+8tt7 2 s 10 (si + 5si 5 ) + 8s 10 (s 2 - 2s 10 si 5 - 2s 2 5 ) 
+m 2 (55m* - 2m 2 (9sio + 46s i5 ) 
+4(-3s 2 + s 10 si5 + 7s 2 5 ))))), 
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S|A 2 | 2 

1 

m 2 m 2 (m 2 -m 2 c + w? w - 2s 10 - 2si 5 ) 4 

x(16(m b + m c ) 2 (Amlm 2 w + 12mjjm c m^ — 6m^m^ 

+(m 2 w - 2(s 10 + s 15 )) 2 (m 4 w - 2m^,(s 10 - 2s 15 ) 
-4si 5 (si + si 5 )) - m^(-13m^ + 8m 2 (s 10 + S15) 
+12(s 10 + si 5 ) 2 ) + m\{2m 2 c m 2 w + 7m 4 w 
+4m 2 w (-5s w + s 15 ) + 4(s 2 - 2s 10 s 15 - 3s 2 5 )) 
+2mjm c (m* - 2m 2 w (13s w + 7«i 5 ) + 8(2s 2 + 3s w s 15 
+sj 5 )) - 2m 2 c (4m 6 w + 4(s 10 + Si 5 ) 2 (3si + 4s 15 ) 
-m^(7si + 4s 15 ) - 8m 2 w (s 2 10 + 3s 10 s 15 + 2s 2 5 )) 
+2m b m c ( y 2m 4 m 2 u - - 16si (sio + s 15 ) 2 
+m^,(-8sio + 4si 5 ) + 4m 2 y (7s 2 + 6S10S15 - «? 5 ) 
+m 2 (-m^ + 2m 2 o (si - 5si 5 ) + 8si 5 (si + S15))) 
+2m t 2 (8mX + 2m^ - 4(s 10 - 4s 15 )(s 10 + s 15 ) 2 
+m^(-llsi + 4si 5 ) + 8m 2 (2s 2 - s w s 15 - 3s 2 5 ) 
-2m 2 (5mt + m^Ssio + Hsi 5 ) 
-2(5 S 2 + 12 Sl0 s 15 + 7s 2 5 ))))), 



1 

m & + m l m c — m l + m c m 2 w + mb{—m 2 c + m 2 — 2si ) 
1 

m 2 m 2 (m 2 - m 2 . + - 2s 10 - 2si 5 ) 3 
x(16(m b + m c ) 2 (5mlm c m 2 u + m^(10m 2 m 2 + 5m^ 

-4s 2 - 14m 2 si 5 + 4sioSi 5 + 8s 2 5 ) + m£m c (-17m 2 m 2 

+ 17m* + 6m 2 (sio - 6s 15 ) - 8(3s 2 + S10S15 - 2s 2 5 )) 

+m 2 m c (-29m^m 2 + 5m^, - 2m^(3si + HS15) 

-4m 2 y (7s 2 - 9s w s 15 - 8s 2 5 ) + 8(3s? + 2s 2 10 s 15 
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-3s 10 4> - 2s 3 15 ) + 2m 2 c (9m 4 w - 2(s w + 2s 15 ) 2 

+m 2 w (3s 10 + 4s 15 ))) - 2m 3 b (2Qm 4 m 2 w - 4m 6 w 

-Um 2 w s 15 (2s w + s 15 ) + m^(9si + 14si 5 ) 

+4(-s? + 2s 2 l0 s 15 + 4s 10 s 2 5 + s 3 5 ) + m 2 c (-15m 4 w 

+2s 10 (5si + 6si 5 ) + m 2 (s w + 18si 5 ))) 

+m c (9m^m^ + m 2 w (m 2 w - 2si 5 )(m^ - 4m 2 si + 4s 2 - 4s 2 5 ) 

+m*(-llm* + 28m 2 u (s 10 + Si 5 ) - 4s 10 (s 10 + 2s 15 )) 

+m 2 c (m b w - 8si Si 5 (sio + s 15 ) - 2m 4 UJ (Qs w + 7si 5 ) 

+4m 2 (3s 2 w + 8si si5 + 6s 2 5 ))) + m b (10m^m^ 

+ - 2sio)(m 2 - 2si 5 )(3m^ - 8771* (sio + S15) 

+4(si + si 5 ) 2 ) - m 4 c (llm 4 w + 8^ + 12s 10 si 5 + 8s 2 5 

-2m^(21s 10 + 25s 15 )) - 2m 2 c (m 6 w + 4mts 15 

+2m 2 w (3s 2 10 - 3s w s 15 - 7sj 5 ) - A(s 3 + 2s 2 s 15 - 4>)))))- 

For W + — > B* + + c + c, we introduce parameters: 

SlO =Pl-P, S15 =Pl-p 5 

with pi, p and p 5 being the momenta of VF + , £>* + and c quark, respectively. 

s|A| 2 

1 

(ml + m 2 m c — mj? + m c m? w + m b (—m 2 c + m 2 — 2si )) 2 
1 

m 2 c m 2 w (m 2 h - m 2 + m 2 w - 2s 10 - 2s 15 ) 2 
x(16(m(, + m c ) 2 (9mlm 2 u + A&m\m c m 2 w + m^(74m 2 m 2 u 

+7mt - 2m 2 w (27s 10 + 2s 15 ) + 4(3s 2 + 3s 10 s 15 + s 2 5 )) 

+2mgm c (27m 2 m 2 y + 5m 4 w - 2m 2 w (55s 10 + 2s 15 ) 

+4(6s 2 + 7s 10 si 5 + 2s 2 15 )) + m 4 (68m 4 c m 2 w - m 6 w 

-8s 2 (3sio + 2s 15 ) - 4mt(7si + 2s 15 ) 

+8m 2 (12s 2 + s 10 si 5 + s 2 15 ) + m 2 (-39mt + 60s 2 

+92s 10 s 15 + 28s 2 5 + m 2 (-254s 10 + 44s 15 ))) 
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-2m h m c {lm\m 2 w - m 2 w (m 2 w - 2s w )(ml + 2m 2 w (s w 
-2s 15 ) + 4(44 - s 10 si 5 + 4)) + m^(-21mt 
+10m^(7si + 2s 15 ) - 4(4 - 5si si 5 + 24)) 
+m 2 c (15m 6 w - 8s w (s w - 2s 15 )s 15 - 2ml(35s w + 12si 5 ) 
+4m^ (144, + 24s 10 s 15 + 34))) + 2m 3 b m c (53m 4 c m 2 w 
-5m 6 w ~ 24s? (sio + s 15 ) - 2m* (13sio + 4s 15 ) 
+12m 2 ,(ll4 + 4) + m 2 ( _ 42m 4 _ 28m^(s 10 - 2s 15 ) 
+4(7s? + 6si si5 + 4s? 5 ))) + m 2 b (54m 6 c m 2 w + (m 2 w 
-2s 10 )(m e w + 4ra* (si - s 15 ) + 8si O si 5 (si + s 15 ) 
+4m 2 w (3s 2 w - 3s 10 s 15 + s 2 15 )) - 4m 2 c (7m 6 w 
-Am 4 w (3s w + s 15 ) + m^(-354 + ^8s 10 s 15 - 2s 2 15 ) 

+4si (4 + S 10Sl5 + 4)) + m c(- 23m t 
+m 2 w (-34s 10 + 52s 15 ) + 4(134 ~ 7sioSi5 + 74))) 
-m\(Y3m\m 2 w + m 4 (-23m 4 w + m 2 w (7As 10 + 28s 15 ) 
+4(34 + 3s w s 15 - Si 5 )) - m l,( m i ~ 2m^(si + 2s 15 ) 
-16si (24 + 4) + Am 2 w (4s? + 3si si5 + 4)) 
+m 2 c (llm 6 w -4m 4 w (17s w + 6s 15 ) 
+4m 2 w (214 + 22s 10 s 15 + 44) + 8(4 + 2s 10 s? 5 ))))), 



S|A 2 | 2 

1 

m 2 c m 2 w (m\ -m 2 c + m 2 w - 2s w - 2si 5 ) 4 

x (16(ra& + m c ) 2 (2mlm 2 w + %m\m c m 2 w + m 4 c (m 4 w - 44 

+4m^,( Sl0 - 2s 15 ) + 44) + (m* - 2(s 10 + s 15 )) 2 (m 4 w 

-2m 2 w (s 10 - 2s 15 ) - 4si 5 (si + s i5 )) - 2m 2 c (m\ 1 

+4si (sio + S15) 2 - ™t( s io + 4si 5 ) - 4m 2 (4 ~ 4)) 
+m 4 b (l2m 2 c m 2 w - m 4 w - 4m 2 w (4s w + s 15 ) + 4(34 + 4s w s 15 
+4)) + 2mlm c (12m 2 c m 2 w - 9ml ~ 6m 2 w (hs w + 3s l5 ) 
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+8(4s 2 + 7s 10 s 15 + 3s? 5 )) + 2m\{\lm\ml - m 6 w 
-12si (sio + si 5 ) 2 + m^(-5si + 4s i5 ) + 4m 2 (5s 2 
+4si si 5 - S15) + m 2 c (-16ml ] - 34m 2 w (s w + s 15 ) 
+44(s 10 + si 5 ) 2 )) + 2m b m c (8m^m 2 - m 6 w 
-16si (sio + si 5 ) 2 + m^(-8si + 4s 15 ) 
+4772 2 o (7s 2 + 6S10S15 - s 2 5 ) + 777 2 (-7777^ 
-2m 2 (5s 10 + llsig) + 8(2s 2 + 5s 10 si 5 + 3s 2 5 ))))), 

SRe(A 1 A 2 *) 

1 

mf + m\m c - m\ + m c m 2 w + m b (-m 2 c + - 2si ) 
1 

m\m 2 w {m\ - m 2 + ml - 2s 10 - 2s 15 ) 3 
x(— 16(m fe + m c ) 2 (4m^m 2 + 21m^m c m 2 w + mjj(38m 2 m 2 

-m A w - 2777^(10510 + si 5 ) + 4(3s 2 + 5s 10 Si 5 + 2s 2 5 )) 

+777ft777 c (39777 2 777 2 - 17777^ - 2777 2 (45Si + 8S15) 
+8(7S 2 + 12S 10 S1 5 + 5S 2 5 )) + 777 2 777 c (35777^ - 15m® 

+ 18777^(510 + S15) - 8(5si - 2si 5 )(si + S15) 2 
+4m 2 (17s 2 - s w s 15 - 8s 2 5 ) + 777 2 (-26777t + 52s 2 
+64si si5 + 32s 2 5 - 6777^(15810 + 8si 5 ))) + 2m\{2§m\m 2 w 
-Ami + 3777t(si + 2si 5 ) + 2m 2 .(8s 2 - 3s 2 5 ) 

+4(-3s? _ 4,2^^ + s 3 s) + m 2 ( _ 1Qm 4 + 3gs 2 o 

+68si si5 + 32s 2 5 - m 2 (69si + 22s 15 ))) + m c {m\m 2 w 
+m 4 c {3m 4 w - 4m 2 w si + 4s 2 - 8s 2 5 ) - m 2 (m 2 - 2si 5 )(m* 

-4777^10 + 4S 2 - 4S 2 5 ) + 777 2 (-3777^ + 8Si Sl 5 (SlO + S 15 ) 

+2ml ] (6s w + 5si 5 ) - 4772 2 o (s 2 + 6si si 5 + 6s 2 5 ))) 
+777 6 ( \Am\m 2 w - (m 2 - 2si )(777 2 - 2si 5 )(3m^ 
-8777^(810 + si 5 ) + 4(si + S15) 2 ) - ^c( m t 
+2tt7 2 (13si + 9si 5 ) - 4(6s 2 + si si 5 - 2s 2 5 )) - 2777 2 (5777^ 
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-4m^(3s 10 + 2s 15 ) + m 2 w (-Us 2 w + Us w s 15 + 22s 2 15 ) 
+4(4 + 2^15-4)))))- 
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